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Suppose that (M, g) is a smooth, compact, connected Riemann manifold of dimension m > 1. 
We denote by \dVg\ the volume density on M induced by g. For any u, G C°°(M) we denote by 
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1. Introduction 

Suppose that (M, g) is a smooth, compact, connected Riemann manifold of dimension 
We denote by \dVg\ the volume density on M induced by g. For any u, G C°°(M) we < 
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{u, v)g their inner product, 

{u,v)g:= / u{x)v{x)\dVg{x)\. 
Jm 

The L^-norm of a smooth function u is then 



Let : C°°(M) — C°°(M) denote the scalar Laplacian defined by the metric g. For L > we set 
UL = UL{M,g):= ker(A - A^), d{L) := dimUL- 

Ae[0,L] 

We equip U l with the Gaussian probability measure. 

d'yi{u) := (27r) ~2~e 2~|c?ti|. 

For any u ^ Ul we denote by 'Nl{u) the number of critical points of u. If L is sufficiently lai^ge, 
then J^l{u) is finite with probability 1. We obtain in this fashion a random vaiiable !N"l = y^L,M,g, 
and we denote by E^Nl) its expectation 



Jut. 



'Ul 

In this paper we investigate the behavior of E^Nl) as L — oo. More precisely we will prove the 
following result. 

Theorem 1.1. For any m > 1 there exists a positive constant C = C{m) such that for any compact, 
connected, m-dimensional Riemannian manifold M we have 

E{'Nl,m,9) ^ C{m)dimULiM,g) as L ^ oo. (1.1) 



The constant C{m) can be expressed in terms of certain statistics on the space the space of 
symmetric m x m matrices . We denote c?7^ the Gaussian measure 'on given by 

(27r) 4 j<j- (^12) 

/i^ = 2(2)+^(m + 2)"^-i. 



Then 



C^M= ( r(l + ^) / |detX|d7*W. (1-3) 



m + 4 / 2 



— '•I'm 

A similai- result holds in the case m = 1. In this case M = and J7l is the space of trigonometric 
polynomials of degree < L. One can show (see [33]) 

R 

-^(■^L,5i) ~ Y 5 L — )• OO. 

We can say something about the behavior of C{m) as m — )• oo. 



'We refer to Appendix B for a detailed description of a 3-parameter family Gaussian measures dTa,b,c on S™, that 
includes ^7, as d-y^ = dTz,i,i- 
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Theorem 1.2. 

TTl 

log C(m) ~ log /m ~ log "I as m — t- oo. (1-4) 

The proof of (1.1) is achieved in two steps. First we prove an integral formula^ (2.2) that expresses 
the expected number of critical points of a function in U l as an integral 

E{JiL)= f PL{x)\dVg{x)\. (*) 
J M 

As explained in [1 1], the integral in the right-hand side of the above equality is the the total curvature 
of the immersion given by the evaluation map 

ev : M Hom(t7i,R), p^ewp, (1.5) 

where evp(u) = u{p), Vit € Ul- 

The integrand pl{x) has a probabilistic interpretation. To formulate it let us denote by Hessa;(w, g) 
the Hessian at x of the random function u ^ U l computed using the Levi-Civita connection of the 
metric g. We identify Hessa;(u, g) with a symmetric linear operators TxM — )■ TxM. Then 

Pl{x) = ^-,.|,\ ^|2x -^( I detHesSa;(tt,g)| \du{x) = o). (1.6) 



Above, E{\du{x)\g) denotes the expectation of the random variable Ul 3 u \du{x)\g S M and 
the quantity 

E ^ I det Hessa; {u,g)\ | du{x) = 

is the conditional expectation of the random variable u i— )• | det HesSa; (it, g) \ given that du{x) = 0. A 
result similar- to (1.6) is proved in [8, 14, 15] using the strategy pioneered by Kac, [25] and Rice, [39]. 
We use a different approach based on the double-fibration trick in integral geometry, [3 1 , §9.1.1]. We 
refer to Subsection 2.1 for a comparison between these two approaches. 

The random vector u i— )• du{x) is Gaussian, and the expectation of |du(a;)|^ can be expressed in 
terms of the covariance matrix of this random vector. Using the regression formula (see [4, Prop. 1.2] 
or (A.2) we express this conditional expectation as the unconditional expectation of a new random 
variable | det ^^(a;)!, where Al{x) denotes a random, Gaussian symmetric m x m matrix whose 
covariance takes into account the coiTclations between the Gaussian variables u i— )• Hessxiu, g) and 
u I—)- du[x). 

Next, we reduce reduce the large L asymptotics of the Gaussian random vector du{x) and matrix 
Al{x) to questions concerning the asymptotics of the spectral function of the Laplacian, i.e., the 
Schwartz kernel of the orthogonal projection onto Ul- Fortunately, such questions were addressed 
in the mathematical literature, [7, 13, 40, 44], by refining the wave kernel method of L. Hormander, 
[22]. The clincher is an asymptotic independence result explained in © below. 

We actually prove a bit more. We show that 

_rn C(m)ujjYl 

lim L 2 p^[x) = , uniformly in a; G M, (1.7) 

where Um denotes the volume of the unit ball in M™. Using the classical Weyl estimates (3.2) we see 
that (1.7) implies (1.1). 

The equality (1.7) has an interesting interpretation. We can think of pL{x)\dVg{x)\ as the ex- 
pected number of critical points of a random function in Ul inside an infinitesimal region of volume 



Formulae of this type appear in the literature under different names, Chem-Lashof or Kac-Rice, depending on one's 
mathematical bias. This author usually resides in the Chern-Lashof part of the mathematical world. 
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|(iV^(a;)| around the point x. From this point of view we see that (1.7) states that /or large L we 
expect the critical points of a random function in U l to be approximatively uniformly distributed. 
We are inclined to believe that as L — )• oo the ratio 

var{yiL) 

has a finite limit g(M, g). Such a result would show that is highly concentrated near its mean value 

as L — )• oo. In [33] we prove that this is the case when M = and moreover, q{S^) 0.4518 

For a holomorphic counterpart of such an estimate we refer to [ ]. The Riemannian case is much 
more challenging since we lack good general off-diagonal estimates for the spectral function. 

We obtain the asymptotics of C(m) by relying on a trick used by Y.V. Fyodorov [18] in a re- 
lated context. This reduces the asymptotics of the integral to known asymptotics of the 1 -point 
correlation function in random matrix theory, more precisely, Wigner's semi-circle law. 

Philosophically, the universality result contained in Theorem 1 . 1 is a consequence of a universal 
behavior of the spectral function £ l along the diagonal. Roughly speaking, if we rescale the metric 
g so that in the limit it becomes flatter, and flatter, then the corresponding spectral function begins to 
resemble the spectral function of the Laplacian on the Euclidean space E™. For a precise formulation 
of this universal rescaUng phenomenon we refer to [26, 34]. 

A related problem was considered by M. Douglas, B. Shiffman, S. Zelditch, [14, 15] where they 
investigate the number of critical points of a random holomorphic section of a large power of a 
positive holomorphic line bundle £ over a Kahler manifold X. In these papers the role of our U l 
is played by the space of holomorphic sections H^{X, L^), and the large L asymptotics is replaced 
by large N asymptotics. The large N asymptotics ultimately follow from the refined asymptotics of 
the Szego kernels obtained by S. Zelditch in [t .: ]. These refined asymptotics then lead to a complete 
asymptotic expansion as — oo for the expected number of critical points of a random holomorphic 
section of L'^ . In our case, the asymptotics of the spectral function £l is very sensitive to the 
background metric g and it is unrealistic to seek a complete asymptotic expansion of ECNl) as 
L — oo that is valid for any metric on g. Finally, the papers [14, 15] do not investigate large 
dimension asymptotics of the type described in our Theorem 1.2. 

The proof of Theorem 1 . 1 reveals several additional interesting universal rescaling phenomena. 

O We identify Ul with = Hom(J7L, M) using the L^-metric. We can thus view the evaluation 
map in (1.5) as a map ew : M ^ U l- For large L this map is an embedding, and we denote hy 
the puUback to M via ev of the L^-metric on U l. Equivalently, if (i/'^,) is an orthonormal basis of 
Ul, then 

k 

The equality (3.9) in the proof of Theorem 1.1 shows that the rescaled metric g{L) := L 2 
converges in the C" topology to Kmg, where g is the original metric on M and Km is a certain, 
explicit constant that depends only on m; see (3.5). This was also observed by S. Zelditch, [44, Prop. 
2.3]. A closely related result was proved in [5, Thm.5]. 

To obtain the convergence of g{L) in stronger topologies we would need bounds on the sectional 
curvature of g{L). In Subsection 3.3 we show that these bounds are equivalent to some refined 
asymptotic estimates satisfied by certain linear combinations of fourth order derivatives of the spectral 
function, (3.20). These estimates hold for homogeneous spaces equipped with invariant metrics. 

A related embedding can be constructed in the holomorphic case and S. Zelditch [ ' -] has proved 
that the resulting sequence of suitably rescaled metrics g^ converges C°° to the original Kahler 
metric. The main reason for such a stronger form of convergence is the better behavior of the Szego 
kernels. Such a regular behavior is not to be expected for the spectral function £l. 
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@ Let us denote by Hessp(it, g) the Hessian at p of a smooth function u defined in terms of the Levi- 
Civita connection V^; see (2.19). Fix an orthonormal basis of TpM so we can identify 'Ressp{u, g) 
with a symmetric m x m matrix. The coiTcspondence 

U L ^ u ^ Hessp(w, 5) G Sm 

defines a Gaussian random symmetric matrix. Similarly, the con^espondence 

Ul^u^ du{p) G T*M 

defines a Gaussian random vector. The estimates (3.13) and (3.14) imply that as L — )• co the random 
vectors 

m+4 m + 2 

L 4 ilessp{u, g) and L 4 duijp) 

converge to independent random Gaussian vectors. Ultimately, this is an asymptotic manifestation of 
the surprising independence result [1, Eq. (12.2. 11) ]. 

© The limit as L — 00 of the random Gaussian symmetric matrix L 2 Hessp(w, g) is a random 
symmetric matrix A^o whose entries a-ij are centered Gaussian variables satisfying the correlation 
equalities 

E{aijaki) = 6ij6ke + Sik5ji + 6ii6jk, \/i,j,k,£. (1.8) 

where Cm is a certain, explicit universal constant that depends only on m. (Up to a rescaling, the 
probability density of is given by (1.2)). Observe that the random matrix A^o, is independent of 
the point p, the choice of metric g and even on the manifold Ml] The universal gaussian ensemble 
Um described by (1.8) is not the GOEm (Gaussian Orthogonal Ensemble) typically investigated in 
random matrix theory and given by the correlation equalities 

EiaijUki) = Sij6ke + dikSje + 6ie6jk, \/i,j,k,i. 

However Um is closely related to GOEm- More precisely a random matrix A in the ensemble Um 
can equivalently be described as a sum of random matrices 

A = B + Ylm, 

where B belongs to GOEm and y is a standard scalar gaussian random variable independent of B. 

The present paper is structured as follows. Section 2 contains the formulation and the proof of 
the key integral formula (1.6), including several reformulations in the language of random processes. 
In this section we also present a simple application of this formula to the number of critical points 
of random spherical harmonics of large degree on 5^. This sheds additional light on a recent result 
of Nazarov and Sodin [30] on the number of nodal domains of random spherical harmonics. More 
precisely, the inequality (2.41) shows that the expected number 5„ of zonal domain on S'^ of a random 
harmonic polynomial of large degree n satisfies the upper bound (5„ < 0.29n^. 

Section 3 contains the proof of the asymptotic estimate (1.1) and Section 4 contains the proof of 
the estimate (1.3). 

For the reader's convenience we have included in Appendix A a coordinate-free, brief survey of 
several facts about Gaussian measures and Gaussian processes in a form adapted to the applications in 
this paper. Appendix B contains a detailed description of a 3-parameter family of Gaussian measures 
on the space of real, symmetric m x m matrices. These measures play a central role in the proof 
of (1.1) and we could not find an appropriate reference for the mostly elementary facts discussed in 
this appendix. Appendix C contains the computations of several Gaussian integrals involving random 
2x2 matrices. 
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Notations 

(i) For any random variable ^ we denote by E{^) and respectively var{^) its expectation and 
respectively its variance. 

(ii) §rn denotes the space of symmetric m x m real matrices. 

(iii) For any finite dimensional real vector space V we denote by V'^ its dual, := Hom( V, M). 

(iv) For any Euclidean space V, we denote by S{V) the unit sphere in V centered at the origin 
and by B{V) the unit ball in V centered at the origin. 

(v) We will denote by <t„ the "area" of the round n-dimensional sphere S"" of radius 1, and by 
u!n the "volume" of the unit ball in M". These quantities are uniquely determined by the 
equalities (see [31, Ex. 9.1.11]) 

7r2 nvra /1\ 

-n^l = nc^n = 2^ = f(rTt) ' = ^' 

where T is Euler's Gamma function. 

(vi) If Vq and Vi are two Euclidean spaces of dimensions no,ni < 00 and A : Vq — > Vi is 
a lineal" map, then the Jacobian of A is the nonnegative scalar J{A) defined as the norm of 
the linear map 

PJ'A : A^Vq A^'Vi, k := min(no, ni). 
More concretely, if ng < ni, and {ei, . . . , e„y} is an orthonormal basis of Vq, then 

J(^) = (detG(^))^/^ (J_) 
where G{A) is the no x no Gramm matrix with entries 

Gij = ( Asi, Asj ) . 

If ni > no then 

J{A) = J{A^) = {detG{A^)f\ (J+) 

where A^ denotes the adjoint (transpose) of A. Equivalently, if dVoli S A^'V* denotes the 
metric volume form on Vi, and dVolA denotes the metric volume form on ker A, then J (A) 
is the positive number such that 

dVolo = ±dVoU A A*dYoh . ( J^) 

2. The key integral formula 

2.1. A Chern-Lashof type formula. As we mentioned in the introduction, a key component in the 
proof of Theorem 1.1 is an integral formula that describes the expected number of critical points as an 
integral over the background manifold M. The goal of this section is to state and prove this formula, 
and then give several probabilistic reformulations. 

The literature on random fields contains many formulse of this type, and their proofs follow the 
strategy pioneered by M. Kac and S. Rice, [25, 39]. As mentioned in the introduction, the general 
integral formula [; , Thm. 4.2] or [14, Thm. 4.4] covers the situation of interest to us. We beheve 
that it would greatly benefit a geometrically inclined reader to see an alternate approach that does not 



He suddenly and untimely passed away in June 201 1. I will miss his generosity and expertise. 
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follow the usual probabilistic pattern pioneered by Kac and Rice, but instead relies on the ubiquitous 
double-fibration trick in integral geometry, [2, 20, 31]. As a matter of fact, our main integral formula 
(2.2) contains as special cases the integral formula of Chem-Lashof, [ ] and Milnor, [ ]. We should 
perhaps describe the differences between the Kac-Rice approach and the double-fibration approach. 

Denote by 5]\j the measure on T*M defined by the integration along the zero section. The Kac- 
Rice approach identifies the measure pg \ dVg | on M as an average (in the space of measures on M) of 
the family of measures {pu '■= {du)*6M'-, u G Ul}, where we denoted by {du)*5M the pullback 
in the sense of distributions of the measure 5m via the map du : M ^ T*M. As is well known, 
[23, Thm. 8.4.2], the pullback operation is well defined provided a certain micro-local transversality 
condition is satisfied. For this reason, the implementation of the Kac-Rice strategy requires two 
intermediary steps, namely replacing the measure 6m with a regulaiization 6^^ followed by a passage 
to limit e \0. 

The double-fibration approach is in some sense dual to the Kac-Rice approach. It leads immedi- 
ately to a description of pg \ dVg \ as the pushforward of a certain measure defined on the total space of 
a fibration ir : P ^ M, and the bulk of the proof is devoted to explicitly computing this pushforwad. 

Suppose that M is a smooth, compact, connected manifold without boundary. Set ?n := dim M. 

Definition 2.1. (a) For any nonnegative integer k, any point p G M and any / G C°°{M) we will 
denote by jkif, p) the A;-th jet of / at p. 

(b) Suppose that U C C°°(M) is a linear subspace. If k is nonnegative integer then we say that U is 
k-ample if for any p £ M and any / G C°°(M) there exists u £ U such that 

jk{u,p) = jkif,p)- □ 

Fix a finite dimensional vector space U C C°°(M) and set iV := dimU. We have an evaluation 
map 

ev = ev^ :M^U'^ := Hom(J7, M), p ^ evp, 
where for any p £ M the linear map evp : J7 — )• M is given by 

evp(ii) = u{p), \/u G U. 

For any u G C°°{M) we denote by N(it) the number of critical points of u. In the remainder of this 
section we will assume that U is 1-ample. This implies that the evaluation map ev*^ is an immersion. 
Moreover, as explained in [32, §1.2], the 1-ampleness condition also implies that almost all functions 
It G J7 are Morse functions and thus J^{u) < oo for almost all u £ U. 

We fix an inner product h = {—,—)fiOnU and we denote by | — |/j the resulting Euclidean norm. 
Using the metric h we identify U with its dual and thus we can regai^d the evaluation map as a smooth 
map ev : M — )■ U. We define the expected number of critical points of a function in U to be the 
quantity 

jiiU ,h) := ^— I J^{u)\dAh{u)\= [ J{{u)^^\dVh{u)\, (2.1) 
(^N-i Js(U) Ju (27r)~ 

V ' 

where <t„_i denotes the "ai^ea" of the unit sphere in M", \dAh\ denotes the "area" density on the 
unit sphere S{U), and |(iV/i('u)| denotes the volume density on U determined by the metric h. A 
priori, the expected number of critical points could be infinite, but in any case, it is independent of 
any choice of metric on M. The space U equipped with the Gaussian probability measure dj^^ is a 
probability space. We denote by J^u the random variable U 3 v ^ yi{u) G Z so that 
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where E{—,d'^f^) denotes the expectation computed with respect to the probability measure d'y^^. We 
will refer to the pair ([/, h) as the sample space. 

The main goal of this section is to express 3sf(J7, h) as the integral of an explicit density on M. To 
describe this formula it is convenient to fix a metric g on M. We will express J^{U,h) as an integral 



Pg{p)\dVg{p)\. 

IM 

The function pg does depend on g, but the density pg (p) \ dVg (p) | is independent of g. The concrete 
description of Pg{p) relies on several fundamental objects naturally associated to the triplet {U , h, g). 
For any p G M we set 

Ul := {u G U; du{p) = O}. 

The 1-ampleness assumption on U implies that for any p £ M the subspace Up has codimension m 
in U so that 

dim Up = N — m. 

Denote by dAgfjjo-^ the area density along the unit sphere S{Up) C U^. 

The differential of the evaluation map at p is a linear map Ap : TpM — U. We will refer to 
Ap as the adjunction map and we will denote by Jg{p) = Jg{p, U) its Jacobian. More precisely, if 
(ei, . . . , Cm) is a (7-orthonormal basis of TpM, then 



Jg{pf = ^ei [ApEi^Apej)^ 



l<i,j<m 



Since ev^ is an immersion we have Jg{p) / 0, G M. 

For any p £ M and any u G U^, the Hessian of at p is a well defined symmetric bilinear 
form on TpM that can be identified via the metric g with a symmetric endomoiphism Hessp(it, g) of 
TpM. We denote this symmetric endomoiphism by Hessp(ti, g). 

Theorem 2.2. If {U, h) is a 1-ample sample space on M, then 

J<{U,h) = ^— / -i- ( / \detHesSpiv,g)\\dAs,uo^{v)\] \dVg{p)\ 
o-N-i Jm Jg[p) \Js{ul) ^ ) 

r I ( r ~ \ (2 2) 

= (2^)"^/ TT:;;^ / IdetHesSpKg)! / \dVu{u)\ \dVg{p)\. 



Proof. Denote by U,^j the trivial vector bundle over M with fiber U, U_j^f := {U x M — >■ M). For 
any p £ M vje. denote by Kp the orthogonal complement of U^ in U. 

Lemma 2.3. The subspace Kp coincides with the range of the adjunction map Ap. 

Proof. Indeed, if {'^n)i<n<N is an orthonormal basis of {U , h), then 



n 



and for any vector field X on M we have 

ApXp = ^(X^„)p^'n. 
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Thus, the function u = Un^n G U, Un E M, belongs to Kp if and only if for any vector field X 
on M we have 

n 

□ 

This proves that the collection {Kp) defines a subbundle K of [/j^/ and the adjunction map induces 
an isomorphism of vector bundle A : TM — )• K. We deduce that the collection of spaces {U^)p^M 
also forms a vector subbundle of the trivial bundle IJ_^ and we have an orthogonal direct sum 
decomposition 

Um = U^®K. 

For any section u of U_}^[ we denote by its i/*^ -component. 

The bundle lJ_^i is equipped with a canonical trivial connection D. More precisely, if we regard 
a section of u of C/^j^ as a smooth map u : M ^ U, then for any vector field X on M we define 
Dxu as the smooth function M ^ U obtained by derivating u along X. The shape operator of the 
subbundle K is the bundle morphism H : TM ® -FC — )• defined by the equality 

■B{X,u) := {Dxuf, yXeC^{TM), ix G 

For every p € M, we denote by Hp the induced linear map Hp : TpM (g) Kp — U^. If we denote 
by Grm{U) the Grassmannian of m-dimensional subspaces of U, then we have a Gauss map 

M 9 p ^ g(p) := Kp £ Gr„(J7). 

The shape operator Hp can be viewed as a Unear map 

Hp : TpM ^ Hom(Kp, J7°) = Tk^ Gr„(J7), 

and, as such, it can be identified with the differential of S at p, [3 1 , §9. 1.2]. Any v G Up determines 
a bilinear map 

Sp-v : TpM ^Kp M, Hp • ^(e,^) := {S,p{e,u),v) ^, 

By choosing orthonormal bases (ej) in TpM and (uj) of Kp we can identify this bilinear form 
with an in x m-matrix. This matrix depends on the choices of bases, but the absolute value of its 
determinant is independent of these bases. It is thus an invariant of the pair {Sp,v) that we will 
denote by | det Hp ■ v\. 

Lemma 2.4. 

J^{U,h) = ^—[ if \deiap-v\\dAsajnAv)\\\dVg{p)\. (2.3) 

(^N-l Jm \Js{Ul) " J 

Proof. Consider the incidence variety 

J ■= {{p,v) e M X S{U); dv{p) = 0} = {{x,v) € M X S{Uy, v^SiU^)}. 
We have a natural double "fibration" 

M ^3 S{U), 

where the left/right projections A, p ai^e the canonical projections. The left projection A : J — )• M 
describes 3 as the unit sphere bundle associated to the metric vector bundle U^. In particular; this 
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shows that J is a compact, smooth manifold of dimension {N — 1). For generic v £ S{U) the fiber 
p~^{v) is finite and can be identified with the set of critical points of v : M ^W. We deduce 

^iU,h) = Atttyt / #p-H'^)\dAHiu)\. (2.4) 

area(5(J7)) Js[u) 

Denote by gj the metric on 3 induced by the metric on M x S{U) and by \dVj\ the induced volume 
density. The coarea formula, [hi, §13.4], implies that 

/ #p-Hv)\dAh{v)\= [jp{p,v)\dV3{p,v)\, (2.5) 
JS{U) Jj 

where the nonnegative function Jp is the Jacobian of p defined by the equality 

p*\dAh\ = Jp-\dVj\. 

To compute the integral in the right-hand side of (2.5) we need a more explicit description of the 
geometry of the incidence variety J. 

Fix a local orthonormal frame (ei, . . . , 6^) of TM defined in a neighborhood in M of a given 
point Pq G M. We denote by (e^, . . . , e™) the dual co-frame of T*M. Set 

fi{p) :=Apei{p) €U, i = l,...,m, p G 0. 

More explicitly, [u) is defined by the equality 

{f,{p),v)^ = deMp), yueU. (2.6) 

Fix a neighborhood It C A~^(0) in M x S{U) of the point (pq, vq), and a local orthonormal frame 
ui{p, v),. . . ,UN^i{p, v) over U of the bundle p*TS{U) M x S{U) such that the following 
hold. 

• The vectors ui{p,v), . . . , Um{p, v) are independent of the variable v and form an orthonor- 
mal basis of K^. (E.g., we can obtain such vectors from the vectors fi{p), . . . , fm{p) via 
the Gramm-Schmidt process.) 

• For {p, v) e U, the space TpEx is spanned by the vectors v),. . . ,UN^i{p, v). 
The collection ui{p), . . . , Um{p) is a collection of smooth sections of U_jyj over 0. For any p £ 

and any e G TpM, we obtain the vectors (functions). 

DeUi{p),...,DpUm{x) G U, 

where we recall that D denotes the trivial connection on C/jvf • Observe that 

3nU = {{p,v) eU; Ui{p,v) = 0, Vi = l,...,m}, (2.7) 
where U-i is the function : x J7 — )• M given by 

Ui{p,v) := {ui{p),v)^. 
Thus, the tangent space of J at {p, v) consists of tangent vectors pQ) i) £ TxM © TvS{V) such that 

dUi{p, i)) = 0, Vi = 1, . . . , m. 

We let ujij denote the m-form 

wc/ := dC/i A • • • A da,n G 17™ (U), 

and we denote by \\uju\\ its norm with respect to the product metric on M x S{U). Denote by \dV\ 
the volume density on M x S{U) induced by the product metric. The equality (2.7) implies that 

\dV\ = \uju A dVE I . 
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Hence 



We deduce 



Jp\dV\ 



\^u\ 



\uju A p*dA\. 



Jp{Po,vo) = Jp{pQ,vo)\dV\{ei, . . . ,em,ui,. . . ,un-i] 



\^u\ 



\uju A p*dS\{ei, ...,em,ui,.. .,un-i) 



1 



\^u\\ L 



Nc/(ei,...,e^ 



l{P0:^'0) 



■^uiPO''"o) 



Hence, 



/ #p~\w)\dA,,{v)\ = f pL.\dV,{p,v) 
Js{u) J3 W^uW 



Sublemma 2.5. We have the equality 



1 



\^u\ 



where J\ denotes the Jacobian of the projection A : J — )• M. 
Proof. Along U we have 



\dV\ 



1 



\^u\ 



\uju A dVj 



1. 



while the definition of the Jacobian implies that 

m = l-\dVg^dAg^^J0J. 

Therefore, it suffices to show that along U we have 

\dV\ = \uju A dVg A dAgijjo^\, 

i.e., 

uju A dVg A d^5.(fj0)(ei, . . . , e,„,ui, . . . ,un-i) 
Since dUi{ui^) = 0, VA; > m + 1 we deduce that 

uju A dVg A (i^5(^0)(ei, . . . , e^, "i, ■ ■ ■ ,un~i) = l^^ui^i, ■ ■ -^Um)]- 
Thus, it suffices to show that 

\uju{ui, ■ ■ ■ ,Um)\ = 1. 

This follows from the elementary identities 

dUi{uj) = {ui,Uj)h = Sij, VI <i,j< m, 
where 6ij is the Kronecker symbol. 

Using (2.9) in (2.8) and the coarea formula we deduce 

/ #p~\w)\dAh{v)\= f if ^u{p.v)\dAsi^^.){v)\] \dVg{p)\. 
Js(u) Jm \ Jsiul) ^ J 



(2.8) 



(2.9) 



□ 



(2.10) 



Observe that at a point {p,v) e X ^ (0) C J we have 

dUi{ej) = {DejUi{p),v)^. 
We can rewrite this in terms of the shape operator Hp : TpM (E> Kp — > Up. More precisely, 

dUi{ej) = {Sp{ej,Ui),v)^. 
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Hence, 

Au{x,v) = \det{Sp{ej,Ui),v)^ \ , 

We conclude that 

#p-\v)\dAh{v)\= I If \detSp-v\\dAs^ijo.{v)\]\dVM{p)\. 
'Sh{U) Jm \Js(ul) J 

This proves (2.3) □ 

To proceed further observe that the left-hand side of (2.3) is plainly independent of the metric g on 
M. This raises the hope that if we judiciously choose the metric on M, then we can obtain a more 
manageable expression for fi{M, V). One choice presents itself. 

Let a be the pullback to M of the metric on V via the immersion ev : M — )■ More concretely, 
for any p ^ M and any X, y E TpM, we have 

ap{X,Y) = {ApX,ApY)^^. 

When g = cr, the equality (2.3) is precisely the main theorem of Chern and Lashof, [' ' ]. 

Fix p G M and a cr-orthonormal frame(ej)i<j<m of TM defined in a neighborhood of p. Then 
the collection Uj = Aej, 1 < j, is a local orthonormal frame of K\q. The shape operator has the 
simple description 

Sp{ei,Uj) = (De^Aej)'^. 
Fix an orthonormal basis (^'n)i<n<Ar of U so that every v ^ U has a decomposition 

a 

Then 

n n 

and 

a 

If G J7p, then the Hessian of at p is a well-defined, symmetric bilinear form Hessp(t') on TpM 
that can be identified via the metric a with a symmetric linear operator 

Hessp(^, a) : TpM TpM. 

If we fix a <T-orthonormal frame (ej) of TpM, then the operator HesSp(t),<T) is described by the 
symmetric m x m matrix with entries d'^^^.v{x). We deduce that 

IdetHp • v\ = \detBesSp{v,a) \ , \/v e S{Up). 
In particular, we deduce that 

y^{U,h) = ^—l I I \detIiesSp{v,a)\\dAg,jjo-.{v)]\dVcr{p)\. (2.11) 
o-N^i Jm \ Js{u%) " J 

Finally, we want to express (2.11) entirely in terms of the adjunction map A. For any p ^ M and any 
V G Up, we define the density 

/ip,^ : K"^TpM R, 



Hp,v{Xi A ■ ■ ■ A Xm) = det{dj^^x,v{p))^^ij^^^ ■ {det{{ApXi,ApXj)h)^^-.^^) 



-1/2 



det(HesSp(^)(X„X,))^^. .^^ • ( det( X,- 



l<i j<m / 
-1/2 



l<i J<»n / 
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for any basis Xi, . . . , Xm of TpM. Observe that for any cr-orthonormal frame ei, . . . , of TpM 
we have 

^p,i;(ei A • • • A Sm) = I det Hessp(i>, a) \ . 
If we integrate fip^^ over v G S{Up), we obtain a density 

\dtiu{p)\ :A"'TpM^R, 

\dnuip)\iXiA---AX^)= / fip^^{XiA--- AX^)\dAs,jjoJv)\, 

Js(Up) 

yXi,...,X.m£TpM. 

Clearly \dfiu{p)\ varies smoothly with p, and thus it defines a density \dfiu{—) \ on M. We want 
to emphasize that this density depends on the metric on U but it is independent of any metric on M. 
We will refer to it as the density ofU. By construction 

J^(U,h) = -^ [ \d^iu{p)\. 

If we now return to our original metric g on M, then we can express \dfiu{—)\ as a product 

\dflu{p)\=Sg{p)-\dVg{p)\, 

where 6g = 5g^u : M — M is a smooth nonnegative function. 

To find a more useful description of pg, we choose local coordinates (x^, . . . , x"*) near p such that 
{d,j.i) is a g-orthononnal basis of TpM. Then 

^^pA^x^ A ••• Aa^„) = det(9^^^.^-j;(p))^<.^^.<^ • [det[{Apdx,,Apdx^)h)^^i j^^)~^''^- 
Observe that the matrix ( d'^.,^.v{p) ) . ^.^^ describes the Hessian operator 

}lessp{v,g) : TpM TpM 
induced by the Hessian of i> at p and the metric g. 

1 /2 

The scalar ( det ( {ApOx^ , Apdxj )h ) ^< j ) is precisely the Jacobian Jg {p) of the adjunction 
map Ap : TpM — > U defined in terms of the metric g on TpM and the metric h on U . We set 

Aa;(V^,5) := / \det^esSx{v,g)\\dAg,jjo-s{v)\. 
Js{ul) " 

Since \dVg{p)\{dx^ A ■ ■ ■ A dx^) = 1, we deduce 

5gy{p)=Ap{V,g)-Jg{p)-\ (2.12) 

This proves the first equality in (2.2). The second equality follows from the first by invoking (A.6) 
and the explicit formula {a) for o-^v-i. □ 

Remark 2.6 (A Gauss-Bonnet type formula). With a little care, the above arguments lead to a Gauss- 
Bonnet type theorem. More precisely, if we assume that M is oriented, then, under appropriate 
orientation conventions, the Morse inequalities imply that the degree of the map p : J — >■ Sh{U) is 
equal to the Euler characteristic of M. If instead of working with densities, we work with forms, then 
we conclude that 

X(M) = (2^)-f / -1^ f / detHesSp(n,g) ^ l_^ \dVh{u)\ \ \dVg{p)\. (2.13) 
Jm ■JgKP) \Jul (2vr)^- / 
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When M is a submanifold of the Euclidean space U , and g = cr, then the above argument yields the 
Gauss-Bonnet theorem for submanifolds of a Euclidean space. □ 

2.2. A formula for variance. The equality (2.2) can be used in some instances to compute the 
variance of the number of critical points of a random function in U . More precisely, to a sample 
space ([/, h) C C°°(Af) we associate a sample space {U^, /ia) C C°°(M x M), the image of U 
via the diagonal injection 

!D : C°°(M) ^ C°°(M X M), u^u^, 
u^{x,y) = u{x) + u{y), \/u£U, x,y £ M. 

The metric h/\ is defined by requiring that the map D : (U, h) — (C/^, /ia) is an isometry. 

There is a slight problem. The sample space is never 1-ample, no matter how large we choose 
U . More precisely, the ampleness is always violated along the diagonal A^/ C M x M. We denote 
by the space of restrictions to M"^ := M x M \ Anj of the functions in U^. Note that for any 
u £ U we have 

y{jjA{u) = ^('u^Ia/i) = X(u)2 - 
If U is sufficiently large, then the sample space is 1-ample and we deduce that 

E{J{jj-y{u,d-fh) =E{J^^jA,d-fi,J. (M2) 

The quantity in the left-hand side of (iW2) is the so called second combinatorial momentum of the 
random variable Jsfjy. Using (2.2) we can express the right-hand side of {M2) as an integral of a 
density over Af^. Since fails to be 1-ample along the diagonal so that the corresponding Jacobian 
term vanishes along the diagonal. We deduce that 

2 

li 




A 6 2 

detHesSp,q(w )| N-2m \dVh{u)\ | \dVg,^g{p,q))\. 
(2-k) 

where U^ ^ ■.= Uln U% 

detB.esSp^q{u^ , g x g) = det Hessp(it, (7) det Hessq(ti, (/), \/u G UpCi Uq, 
and is the Jacobian of the adjunction map ^ : Tf^^^q^M x M — ;> U^. 

2.3. A Gaussian random field perspective. The formula (2.2) looks hopeless for two immediately 
visible reasons. 

• The Jacobian Jg (p) seems difficult to compute. 

• The integral Ip in (2.2) may be difficult to compute since the domain of integration Up may 
be impossible to pin down. 

We will deal with these difficulties simultaneously by relying on some probabilistic principles in- 
spired from [ i ]. For the reader's convenience we have gathered in Appendix A the basic probabilistic 
notions and facts needed in the sequel. 

Consider again the metric <t = au, the puUback of the metric honU via the evaluation map. We 
wiU refer to it as the stochastic metric associated to the sample space {U,h). It is convenient to have 
a local description of the stochastic metric. 

Fix an orthonormal basis Vi, • • • , iPn of The evaluation map ev*^ : M — )• Z7 is then given by 

M3x^Y,'il^^{x)-'il^^€U. 
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If p G M and U is an open coordinate neighborhood of p with coordinates x = (x^, . . . , x™), then 
M5.m9..)=E^(p)^(p)' Vl<.,j<m. (2.14) 

n 

Note that if the collection (9j.i)i<i<m forms a (j(-orthonormal frame of TpM, then 

>/g(p)^ = det ap{d^i,d^j) . (2.15) 

L J l<i,j<m 

To the sample space {U, h) we associate in a tautological fashion a Gaussian random field on M 
as follows. The measure d'^i^ in (2.1) is a probability measure and thus {U,d'yh) is naturally a 
probability space. We have a natural map 

C : M X U ^ R, M X U B {p,u) ^ ^p{u) := u{p). 

The collection of random variables {Cp)peM is a Gaussian random field on M. 
Using the orthonormal basis (i/'^) of U we obtain a linear isometry 



3t = {ti, . . . ,tn) ^ Ut = "^tklpi, G U, 



(2.16) 



with inverse u i— )• ifc('u) = h{u, ■0^.). For any p ^ M and any t G we have 

k 

The covariance kernel of this field is the function £ = £[7:MxM— ^-M given by 

£(p, q) = E{^p, = V f / t,tfcfi7iv(*)) tPj{p)tPk{q) 
M 
A;=l 

where d7Ar is the canonical Gaussian measure on M^. 

If p G M and U is an open coordinate neighborhood of p with coordinates x = {x^, . . . , x"*) such 
that x{p) = 0, then we can rewrite (2. 14) in terms of the covariance kernel alone 

5^£(x, y) 

0-p(5x«,9^j) = -^^-^^U=y=0- (2.17) 

Note that any vector field X determines a new Gaussian random field on M, the derivative of u along 
X. We obtain the Gaussian random variables 

u I— )• {Xu)p, u I— (Yu)p, 

and we have 

cTp{X,Y) = E{{Xu)p,{Yu)p). (2.18) 

The last equality justifies the attribute stochastic attached to the metric a. 

We denote by V the Levi-Civita connection of the metric g. The Hessian of a smooth function 
/ : A/ — M with respect to the metric g is the symmetric (0, 2)-tensor V^/ on M defined by the 
equality 

V^f{X,Y) := XYf - {VxY)f, yX,Y G Vect(A/). (2.19) 
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If p is a critical point of / then Vp/ is the usual Hessian of / at p. More generally, if (x^, . . . , x*") 
are g-normal coordinates at p, then 

V^/(5,.,5,,) = dl^J{p), VI < i,j < m. 

For any p S M and any / G C°°(M) we use the metric gp to identify the bilinear form V^/ on 
TpM with an element of §>{TpM), the vector space of symmetric endomorphisms of the Euclidean 
space {TpM, gp). For any p G Af we have two random Gaussian vectors 

U Bu^ Vlu G S(TpM), UBu^ du{p) G T*M. 

Note that the expectation of both random vectors are trivial while (2.17) shows that the covariance 
form of du{p) is the metric a p. 

To proceed further we need to make an additional assumption on the sample space U. Namely, in 
the remainder of this section we will assume that it is 2-ample. In this case the map 

U Vlu G §>iTpM) 

is surjective so the Gaussian random vector VpW is nondegenerate. A simple application of the co- 
area formula shows that the integral Ip in (2.2) can be expressed as a conditional expectation 

Ip = ^^ddetV^wl \du{p) = 0). 

Observing that 

Jg{p) = {detSauip))K (2.20) 

we deduce that 

^{U,h) = -^ [ {detSduip))-'^E{\detVlu\\du{p) = 0)\dVg{p)\. (2.21) 

The last equality is very similar to the main conclusion of the Expectation Metatheorem, [ 1 , Thm. 
11.2.1] or the expectation formula in [4, Thm. 6.2]. We can simplify the equality (2.21) even more by 
taking full advantage of the Gaussian nature of the various random variables involved in this equality. 
The covariance form of the pair of random variables VpU and du{p) is the bilinear map 

n : ${TpM f X TpM M, 

n{^,v) = E{{^,Vlu)-{du,r])), y^eSl, 11 TpM. 

Using the natural inner products on ${TpM) and TpM defined by gp we can regard the covariance 
form as a linear operator 

rtp : TpM §>{TpM). 

Similarly, we can identify the covariance forms of VpU and du with symmetric positive definite 
operators 

5v2„ : S(TpM) ^ S(TpM) 

and respectively 

Sduip) ■■ TpM ^ TpM. 
Using the regression formula (A.4) we deduce that 

£;(|detV^it| I du{p) = O) = E{\d(iiYp\), (2.22) 

where Yp -.U ^ §{TpM) is a Gaussian random vector with mean value zero and covariance operator 

Sp = Sy, := Sv2^ - f^S-J(p)f^t . g(^rpM) ^ SiTpM). (2.23) 
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Since U is 2-ampIe the operator Hp is invertible and we have 

dimS(Tp) 1 f (Sp^Y.Y) 

£;(|detyp|) = (27r) 2 — {detSp)-2 |dety|e 2 dVg(Y). (2.24) 

Js(Tr,M) 



l§{TpM) 

We deduce that when U is 2-ample we have 

J{{U,h) = —^f (detSrf,(p))^^£;(|detyp|)|dF,(p)|, (2.25) 
(27r) 2 Jm 

where is a Gaussian random symmetric endomorphism of TpM with expectation and covariance 
operator Hp described by (2.23). 

To compute the above integral we choose normal coordinates {x^ , . . . , x") near p and thus we can 
orthogonally identify TpM with M™. We can view the random variable VpU as a random variable 

i/P : (7^S„ :=S(IR'"), U 3 u ^ HP{u) Hf.{u) = d^^Mp), 

and the random variable du{p) as a random variable 

DP -.U ^ M'", u ^ DPu G R"", Dfu = d^.u{p). 

The covaiiance operator Sdu{p) of the random variable is given by the symmetric m x m matrix 
with entries 

(9^£ (x, y) 

ap{d^^,d^,) = \^=y=Q. (2.26) 

To compute the covariance form S//p of the random matrix we observe first that we have a 
canonical basis {£,ij)i<i<j<m of so that ^ij associates to a symmetric matrix A the entry aij 
located in the position Then 



N 



(2.27) 



_^ dx^dxi dy^dy^ a: j/ 

Similarly we have 



= E{dl.^Mp),d,.u{p)) = £|^U=.=o- (2-28) 

To identify Q with an operator it suffices to observe that (d^k ) is an orthonormal basis of TpM, while 
the collection { }i<j C S^, 

\V2^i„ i<j 

is an orthonormal basis of S^. If we denote by Eij the dual orthonormal basis of S^, then 

^d^^ = X]*^(^'*J'^^'■■)^^J■ 
Remark 2.7. If the metric g coincides with the stochastic metric a, then the covariance operator U is 
trivial. For a proof of this and of many other nice properties of the metric a we refer to [1, §12.2]. □ 
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2.4. Zonal domains of spherical harmonics of large degree. In the conclusion of this section we 
want to discuss an immediate application of the above results to critical sets of random spherical 
harmonics. 

Let (M, g) be the unit round sphere 5^. The spectrum of the Laplacian on is 

A„ = n(n + 1), n = 0,1,2,..., dimker(A„ — A) = 2n + 1 = d„. 

The space Un = ker(An — A) has a well known descrition: it consists of sperical harmonics, i.e., 
restrictions to 5^ of harmonic polynomials of degree n in three variables. We want to describe the 
behavior of 3sf([/„) as n — )• oo, where J7„ is equipped with the L^-metric. In other words we want 
to find the expected number of critical points of a spherical harmonic of very large degree. 

In this case the covariance kernel £n(p, q) of Un has a very simple description. More precisely, if 
(^'fe)i<fc<2n+i is an orthonormal basis of C/„, then the classical addition theorem, [29, §1.2] shows 
that 

k 

where • denotes the inner product in M^, and P„ denotes the n-th Legendre polynomial, 

1 r/" 

= (-1)"^ (1 _ t^)^. 
^ ' ^ ' 2"n!dt"' ^ 

In this case the stochastic metric cr = <t„ is obviously S'0(3)-invariant and it is a (constant) multiple 

of the round metric. In view of Remark 2.7 this implies that for any p £ the random variables 

Un B u Hessp(ii,5r) and Un B du{p) 

are independent and we deduce that 

/ \ 

Js2 Jg{p) Ju„ 27r^— 



\ =:rf7n('") / 



Clearly, the integrand in the above formula is invariant with respect to the 50(3)-action on and 
we thus have 

^iUn) = ^^ [ \det-ResSp^{u,g)\d^n{u), (2.29) 
JgKPo) JU„ 

where pg a fixed (but arbitrary) point on 5^. To compute the term in the right-hand side of the above 
equality we use the equalities (2.26) and (2.27). 

Fix normal coordinates {x^,x'^) in a neighborhood of Pq so we can view £„ as a function 
E.n{x, y). The location of a point p G is described by a smooth function 

B {x^,x^)^p{x^,x^) G M^. 

The tangent vector d,j.i, viewed as a vector in M?', corresponds with the derivative p^i := d^ip of the 
above function. At p^ we have 

Px' 'PxJ = ^ij and p^., .po = 0, (2.30) 

The arcs Ci = {x^ = 0} and C2 = {x^ = 0} are portions of great circles intersecting orthogonally 
at Pq. Note that is the arclength parameter along Ci, i = 1,2. The vectors p^i ai"e unit tangent 
vectors along these arcs. This shows that at Pq we have 
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Since the arcs Ci and C2 ai^e planar their torsion is trivial and the Frenet formulae imply that at Pq we 
have 

Pxixi =0, Vi / j. 

The last two equalities can be rewritten in compact form as 

Px^xo = -SijPo, (2.31) 

We set 



2"'+lw/,N 2n + l n{n + l) 1 , 

._ 2n+l (2n + l) (n + 2)(n + l)n(n-l) _ 1 , 

" ■ TT " ^ 47r 16 327r ■ 



(2.32) 



p=qr 



( Pn^ {P'Q) {Px- • V ) (P^;^ • ) + ^i^^ (P • Q) (Px^ • Qy« ) (Px« • Qj/O ) 



We deduce 

= ~~r^ ( (P • Q)Pxi • Qyfe + Pn^ (P • 9) (Px-J • Q){P • Qy'' )) (2-33) 

and 

-^9(Po) = (2.34) 
To compute 5^^^^^^^^ £n(p, q) at p = q = pg we will use (2.30) and (2.31) to cut down the 
complexity of the final formula. We deduce that at p = g = we have 

dli^jykyl^n{P, q) = '^^^ ^ [P'niP • q)Px'x^- • Vy'-' + ^nHP • Q){Px^x1 » Q)iP » Qy^'y'^) 

2n + l 
An 
and thus 

d^ixiykye^n{P,q)p=q = {Sn +tn)SijSke + tn{SiiSjk + SikSji). (2.35) 

Denote by dVn the pushforward of the Gaussian measure dj^ via the Hessian map 

Un3 Hesspg (u, 5) G S(rppS'^) = 82- 
We deduce from (2.35) that the covariance form I]„ of dVn satisfies the equality 

where Sa,b,c is defined by the conditions (B.2a) and (B.2b). Observe that a„, c„ satisfy (B.4), 

— ~t~ 2C72 . 

As explained in Appendix B, this implies that dTn is O (2) -invariant. Set 

* (^n , * * 

ijn ijn ijn 

and denote by dF* the Gaussian measure on §2 with covariance matrix Sa*,b* ,c* ■ Using (A.7) we 
deduce that 

/ \detX\dTniX) = tn I \detX\dTl{X). 
From (2.29) and (2.34) we now deduce 

2t I' 

X(J7„) = — / |detX|dr;(X). (2.36) 
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Observe that as n — )• oo we have 
SO that 

Ji{Un) ~ — / I detXl dr3,i,i(X), (2.37) 

where dFa is the Gaussian measure on §2 with covariance form 1 1. More precisely (see 

(B.ll)) 



4(27r)l .^h^. 



dx 2 J • 



l<j<J<2 



In Appendix C we show that 

/ |detX|dr3,i,i(X) = (2.38) 
J%2 v3 



and we deduce from (2.29) that 



Ji{Un) 1= as n ^ 00. (2.39) 

v3 



Let us observe that for n very large, a typical spherical harmonic u G Un is a Morse function on 
5^ and is a regular value. The nodal set {u = 0} is disjoint union of smoothly embedded circles. 
We denote by the set of connected components of the complement of the nodal set are called 
the nodal domains of u and we denote 5{u) the cardinality of D^. A result of Pleijel and Peetre, 
[6, 35, 38], shows that 

5(u) < w 0.692n2, (2.40) 

where jo denotes the first positive zero of the Bessel function Jq. 

We think of 6{u) as a random variable and we denote by 5„ its expectation. 



J U n 



" (2^)diml 

Recently, Nazarov and Sodin [30], have proved that there exists a positive constant a > such that 

5n ~ an^ as n — 7- 00. 

Additionally, for large n, with high probability, 5{u) is close to a-n? (see [30] for a precise statement). 
The equality (2.40) implies that 

4 

a < — w 0.692. 
jo 

We can improve this a little bit. 

Denote by p{u) the number of local minima and maxima of u, and by s{u) the number of saddle 
points. Then 

yi{y) = p{u) + s{u), p{u) - s{u) = x{S^) = 2. 

This proves that 
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For every nodal region D, we denote by p{u, D) the number of local minima and maxima"* of u on 
D. Note that p{u, D) > for any D and thus the number p{u) = J2De'D ^) viewed as 

a weighted count of nodal domains. We set 

PiUn) := ^ Lu. [ e-^l-I^H \du\. 
(27r) — 2 — Jc/„ 

The equahty (2.39) imphes that 

p{Un) ~ — i=n'^ as n — )• oo. 

This shows that 

a < ^ 0.288. (2.41) 

2\/3 

Remark 2.8. Using (2.13) we deduce in a similar fashion a stochastic Gauss-Bonnet formula 

2 = x('S'^) = — I detHesSp,j(it, d-fn{u). 

Sn Ju„ 

In Appendix C we give a direct proof of this equality to test the correctness of the various normaliza- 
tion constants in the above computations. □ 

3. The PROOF OF Theorem 1.1 

3.1. Asymptotic estimates of the spectral function. We fix an orthonormal basis of L^(M, g) con- 
sisting of eigenf unctions of A^, 

Ag^-n = Xn'^n, 71 = 0, 1, ... , Aq < Al < • • • < An < • • • . 

The collection (^n)A„<L is therefore an orthonormal basis of Ul so that the covariance kernel of the 
Gaussian field determined by C/^ is 

£l(P,9)= J2 ^n(p)^n(Q). 

A„<L 

This function is also known as the spectral function associated to the Laplacian. Equivalently, £ l can 
be identified with the Schwartz kernel of the orthogonal projection onto Ul- Observe that 



/ £L{p,p)\dVg{p)\=dimUL. 

J M 



In the groundbreaking work [22], L. Hormander used the kernel of the wave group e**^ to produce 
refined asymptotic estimates for the spectral function. More precisely he showed (see [22] or [24, 
§17.5]) 

£l(p,p) = (^^^ + OiL"^ ) as L ^ oo, (3.1) 

uniformly with respect to p G M. Above, denotes the volume of the unit ball in M"\ This implies 
immediately the classical Weyl estimates 

d\mUL ~ ^^vol3(M)Lf . (3.2) 



simple application of the maximum principle shows that on each nodal domain, all the local extrema of y are of the 
same type: either all local minima or all local maxima. Thus p{u, D) can be visualized as the number of peaks of |u| on 
D. 
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Hoimander's approach can be refined to produce asymptotic estimates for the behavior of the deriva- 
tives the spectral function in a neighborhood of the diagonal. We describe below these estimates 
following closely the presentation in [7]. For more general results we refer to [40, Thm. 1.8.5, 1.8.7]. 

We set A := La . Fix a point p and normal coordinates x = (x^ , . . . , x"^) at p. Note that x{p) = 0. 
For any multi-indices a, /3 G Z>q we have (see [7, Thm. 1.1, Prop. 2.3]) 

d°'+^EL{x,y) . , . yn+\a\ + \l3\ y^m+\a\+\l3\-l\ 

dx^dyl^ U=y=o - <^mia,PM -\-uyA ), 

where 

' a - /3 (2Z)'" 



Cm(a,/3) = < 



/b™ x^+^\dx\, a-f3£ (2Z)" 



(277)" 

and B"^ denotes the unit baU 

B"^ = {xe M™; \x\ = 1}. 
The estimates (3.3) are uniform in p G M. Using (A.6) we deduce (compare with (B.13) ) 

X '^'^\dx\ = ^m_, i„i4-ifli4-^N / ^ '^^ —\dx\. 



We set 

Km = Cm{a,a), \a\ = 1 

so that 

1 



vr 2 



(3.3) 



(3.4) 



^"^ " (4^)f r(2 + f ) L vrf " 2(4vr)f r(2 + f ) " ^^"^^ 
For any i < j define aij G so that 

a;"*^' = XiXj. 

For i < J and k < i we set 

1 /■ e"l^l^ 

Crn{iJ]kJ)=Crn{aij,ake) = ——^— — / XiXjXkXi — j^\dx\. (3.6) 

(47r) 2 r(3 + Jr™ 



For i < j we have 

1 /" 2 2^~'^'^ 1 

Cmihi] jj j) = 7- ^ iri_ /„ ^;7T / ^i^i = — — m , — r- =: Cm- (3.7) 
(47r) 2 r(3 + f ) Jr™ 7r 2 4(47r) 2 r(3 + f ) 

Finally 

1 [ 4e~l^l^ 3 

Cmih^h^) — ~~. , ITT / I^l'^'^l ~ TT^i X ™ /o mT ~ (3-8) 

(47r) 2 r(3 + f ) Jr- 7r 2 4(47r) 2 r(3 + f ) 

and 

C^(i,j;A:,^) =0, VA: < ^, (i, j) ^ (A:,^). 
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3.2. Probabilistic consequences of the previous estimates. We denote by cr^ the stochastic metric 
on M determiner by the sample space Ul, L ^ 0. As explained in Subsection 2.3 the covariance 
form of the random vector Ul 3 u du{p) G T*M is <Tp, and from (3.3) we deduce 

= K,nX"'+^9p{d^^,d^j) + 0{X"'+^) as L ^ oo, uniformly in p. 

In particular, if -S^„(p) denotes the covariance operator of the random vector du{p), then we deduce 
from the above equality that 

Sduip) = K^y^-^^lm + 0(A™+i), uniformly in p, (3.10) 
and invoking (2.20) we deduce 

r r 1 — m(m + 2) , m(m + 2) ,^ 

J^{p) = (det 5i(p)) 2 = K;h + 0{\— 1) , uniformly in p. (3.11) 

Denote by S|^p the covariance form of the random matrix 

Using (2.27) and (3.3) we deduce 

I]^p = c„A™+^S3,i,i + 0(A'"+3), uniformly in p, (3.12) 

where the positive definite, symmetric bilinear form "^3,1,1 : x — )• M is described by the equal- 
ities (B.2a) and (B.2b). We denote by T^ i i the centered Gaussian measure on S,„ with covariance 
form S3,i,i- 

The equahty (2.28) coupled with (3.3) imply that the covariance operator ftp satisfies 

ft^ = 0(A"*+2), uniformly in p. (3.13) 

Using (3. 10), (3.12) and (3.13) we deduce that the covariance operator Hp defined as in (2.23) satisfies 
the estimate 

= c™A™+^Q3,i,i + 0(A"*+^), as L^oo, uniformly in p, (3.14) 

where ^3,1,1 is the covariance operator associated to the covariance form Il3,i,i and it is described 
explicitly in (B.3). If we denote by cITl the Gaussian measure on S,„ with covariance operator H^, 
we deduce that 

(2^) — (det 3^) 2 f<^ 

' V ' 

\dY\ 

where 

AT J- c m{m + l) 

Nm. = dim brr, = 



2 

Let us observe that \dY\ is the Euclidean volume element on defined by the natural inner product 

on §rn, {X, Y) = tr{XY). We set 



cl • — Cm A , Q„ — '— '!)• 

^ CL 



Using (A.7) we deduce that 



1 f (ct)'- 

' 'dety|e —\dY\ = -I \detY\e "2 — \dY\. 



(27r)^(detH^)2 Js„, (27r)^ (det Q^)3 
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From the estimate (3.14) we deduce that 



1,1 as L -)• oo, uniformly in p. 



We conclude that 



r EL m(m+4) f 

E{\detYp\)= \ det Y\dT LiY) c^X^^ | dety|dr3,i,i(y). (3.15) 
The measure dFs 1 1 is described explicidy in (B.l 1), more precisely 

<ir3,i,i(y) = ^ e-i(-^^-^(-^)0|dy|, 

(27r) 2 ^j-fi;;; 

where /x^ is given by (B.12). Using (2.25), (3.11) and (3.15) we deduce that 

\ 2 m(m + 4) »n(m + 2) (' 

j^) x^-^ — ^voi,(M) / |dety|(ir3,i,i(y) 

(3.2, /|j,\9(2^r,i,„^^ 

Observe that 



Cm_ _ r(2+f) _ _ ttT (27rr _ ,^^(, ^ m 

i^™-2r(3 + f)-m + 4' ^'"-^l + f) um -^4^^^H'+2 

This completes the proof of (1. 1) and (1 .7). □ 

3.3. On the asymptotic behavior of the stochastic metric. We denote hy g{L) the metric 

g(L) := A-(™+2)^L ^ L-'^^L^ 
where Km is described by (3. 5). The estimate (3.9) shows that 

g{L) — > g as L oo, 

where Km is described by (3.5). The metrics g{L) are closely related to the metrics constructed in 
[5, Thm. 5]. We want to discuss here possible ways to improve the topology of the convergence. 

Observe that if g{L) were to converge in the C^-topology to Km then the sectional curvatures 
of g{L) would have to be uniformly bounded. Conversely, the results of S. Peters [16] show that 
the convergence coupled with an uniform bound on the sectional curvatures would yield a C^'" 
convergence. 

The results in [ ' , §12.2.1] describe a simple way of expressing the sectional curvatures of in 
terms of the spectral function El- Here are the details. 

Denote by the Levi-Civita connection of the metric a^. Fix a point p ^ M and g-normal 
coordinates {x^, . . . , x"^) at p. We set 

h,...,ia;ji,..;jb ■ Qrj-n . . . dx'-'^dyi^ ■ ■ ■ dy^b 



cr{L)ij := ap{d^r,d^j), 1 <i,j < ?n, 
„ the inverse matrix of ( <T(L)jj )i< 



and we denote by {a{Ly^ )i<i,j<m the inverse matrix of {(j{L)ij )i<i,j<m- From [I, Eq. (12.2.6)] 
we deduce 
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We set 

e I 

so that 

k 

For u G J7l we set 

Hl^iu) := {d9.d,,u - iVi)J,,)u)^ = ds.d,Mp) -YnL%d.,,u{p)- (3-16) 
We think of the matrix Hf-iu) as an element H^{u) ^T*M® T*M, 

and we set 

H^{u) A H^{u) := Y Htjiu)H^iiu)dx' A dx'' dx^ A dx^ 
= : Y Qikjli:"')dx' A dx^ (g) dx^ A dx^. 

i<k,j<t 

Note that 

Qi^,{u) = 2iHt,iu)Hteiu) - Hj;^{u)Hk{u)). 
We denote by the Riemann tensor of cr^ and we set 

The map 3 u Qikjii'^) £ 1^ is a random variable and according to [1 , Lemma 12.2.1] we 

2i?4^., = -£;(gf,^.,). (3.17) 



have^ 



In particular we deduce that 

From (3.9) we deduce that 

c7{L)ij = ef^j ~ Kr.^X^+'^Sij + 0(A"'+^) as L ^ oo. 

Hence 

From (3.3) we deduce that as A — cxo we have 

- E -k;^2 i^"' + 0(A-^) ) £^-, - K;;^^'r,k + 0{X'') = 0(1), (3.18a) 

E{dl.^Mp),d,^u{p) ) = 8^^, = 0(A-+2) (3.18b) 
Using the estimates (2.27), (3.16), (3.18a) and (3.18b) in (3.17) we deduce 

^{^iiHj^j - {H^f ) = {£-ii-jj - S-ij-ij ) + 0(A'"+^). 



^Alternatively, in our case, the equalities (3.17) are simple consequences of Theorema Egregium, [31, §4.2.4, Eq. 
(4.2.12)]. 
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We deduce that the sectional curvature of along the plane spanned by d^^i , d^k is 

jy-L ^ijij ^ (pL _ pL \ \ o \ ~ ^ir,ii \ 

On the other hand 

where Cm{i, j] k, £) is defined by (3.6), and we deduce 

- = {Cn.{i,rJ,j) - Cm{i,j;i,j))X'^+^ + 0(A™+3) = 0(A'"+3). (3.19) 

Hence 

^ij ^ \2m+4 i^iijj ~ + ^("^ )" 

The sectional curvature of g{L) = A~'"~^cr^ along the plane spanned by d^i,d^j is 

TC^ = \™+2 — (pL _ pL \ , (O/i N 

We deduce that the sectional curvatures of g{L) are uniformly bounded if and only if 

£.fi.jj - efj.ij = 0(A™+2) uniformly over M. (3.20) 

Note that the estimates (3.20) are stronger than the estimates (3.19) which are direct consequences of 
the Bin-H6rmander estimates (3.3). 

Let us point our that (3.20) hold when (M, g) is a homogeneous space equipped with an invariant 
metric. Indeed, in this case the metric g{L) has the same symmetries as g and thus there exists a 
constant cl > such that g{L) = cig. Then — )• Km as L — )• 00 so that g{L) — )• K^g in the 
C°° -topology and therefore K^j = 0(1). 

4. The proof of Theorem 1.2 

4.1. Reduction to the classical Gaussian orthogonal ensemble. We begin by describing the large 
m behavior of the integral 

'...^■=^^ / |detX|.-H"-=^<")')|,X|. 

(27r) 4 ^//^ J^m 

where we recall that 

^^ = 2(2)+™-i(m + 2). 
We will use a trick of Fyodorov [18]; see also [17, §1.5]. Recall first the classical equality 

For any real numbers u,v,w,we have 

ut'^ + V tr{X + wtlmf = (n + mnP')t^ + 2vw{iT X)t + v tr 

=: a{u, V, w)t^ + 6(n, v, w)t + c{u, v, w). 

We seek u,v,w such that 

v'^w'^ . ,,^9 V 6^ — 4ac 1 / 1 



„-(trX)^ ^trX^ = = --{ tiX^ (trX)' 

u + mw^ u + mw^ 4a 4 \ m + 2 



RANDOM SMOOTH FUNCTIONS 



27 



We have 



and we deduce 



Hence 



V 1 v'^w'^ 1 

u + mw'^ 4' u + mw"^ 4(m + 2) ' 

vw^ = -, V = -(u + mt(;^)<^=>u = — muP' . 

(m + 2) 4 



1 , m 



vo^ = — u = Av 



v{m + 2) v{m + 2) 



We choose = ^ so that 



2 2m 4 , , 

u; =7 r, u = 2—- - = , a(u,v,w) = 4v = 2, 

(m + 2)' (m + 2) m + 2' ^ ' ' ^ 



1 



TT 



2 

d7(s) 



Hence 



22 (27r) — 4 — ^//^ -^K V^S^ / 



(4.1) 

/ ( / \det{xl^-Y)\e-"2''^'\dY\] djix). 



= -.fm{x) 

For any 0(n) -invariant function : S„ — )• M we have a Weyl integration formula (see [3, 17, 27]), 

1-^ [ f{X)\dX\ = ^f /(A)|A,„(A)||dA|, 

(27r) 2 JS„ -^n Jr" 

where 

An(A):= n (^^■-^^)' 

l<2<j<n 

and the constant Z„ is defined by the equality [3, Eq. (2.5.11)] , 

„ n 

Z^:= e-|l^l>„(A)||dA| =2tn!r[r(^). (4.2) 
Now observe that for any Aq € M we have (with defined in (4.1)) 

. dim ^rn 2 I ^ \ 

/m(Ao) = ^^""t ' / IT 1^^- - -^Ol I ^-(^) I l^^l 
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1^2 , dim Sm 



= / e 2l^,=i\\Am+i{Xo,Xi,...,\m)\\dXi---dX^\ 

-X^ / \ dim Sm 

= ^'"^'"^^^ ^--^^ / e-l Sr^i I A^^.i(Ao, Ai, . . . , A^) I IdAi . . . dA^| . 

^ V ' 

=:Pm+i(Ao) 

The function Rn{x) = npn{x) is known in random matrix theory as the 1-point correlation function 
of the Gaussian orthogonal ensemble of symmetric n x n matrices, [13, §4.4.1], [19, §3], [27, §4.2]. 
We conclude that 

, . dim Sm . . I 

Im = / Pra+l{.x)e 2 d'y[x) = / Pm+l[,x)\ -C 2 dx. 

^m. . hs ^m. .ha V vr 



We have 



22(m+ l)r 



m + 1 



dim Srn I / ^ , > 1 



22 (27r) 4 



vr \ 2 J J^. 



V 2 y 22I 2 J+^~(m + 2)2 
2, , r(^^^) /2" 2T(^) 



22I 2 V vr22l 2 



We deduce 



We set 



and we deduce 



2 2r ( "^+3 ^ 

2 ^ - / pm+i(x)e~'^dx. (4.3) 



7r22(, 2 )+— 



Pn{s) ■■= y/npniVns), 



Pn{x)e '2 dx = pn{\/ns)e 2 ds = n e 2 pn[s)ds 
Jr Jr 

27r\ 3 f (3n)3e-^ _ , , , (4.4) 

^ / ——1 Pn{s)ds. 

3ny 7k (27r)2 

= :Wn{s) 

To proceed further we use as guide Wigner's theorem, [3, 13, 17, 27] stating that the sequences of 
probability measures 

Pn{x)dx = y/npn{-\/nx)dx = —=Rn{\/nx)dx 

\ n 
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converges weakly to the semi-circle probability measure^ p{x)dx, 

= (4.5, 

^ ' vr |0, \x\ > y/2. 

We observe that the Gaussian measures Wn{s)ds converge to the Dirac delta measure concentrated at 
the origin. This suggests that 

f _ 

lim / pn{s)wn{s)ds = p{0) = — . (4.6) 

n^oo vr 

We will show that this is indeed the case by slightly refining the arguments in one particular proof 
of Wigner's theorem; see [ I 7 , §7. 1.6], [19, §6.1] or [27, A.9]. For the moment we will take (4.6) for 
granted and show that it immediately implies (1.4). 
Using (4.6) in (4.3) and (4.4) we deduce that 

^ 2r(^) / 27r Y V2 

1 An-n rn-i TT7 TT ^ as m OO. 

We now invoke Stirling's formula to conclude that 

f m + 3\ m 

log Im log r ( — ^ — I ~ y log m, as m oo. (4.7) 

Form (1.3) we deduce that 
log C{m) ■ 
Stirling's formula and (4.7) imply that 



log C(m) = log + — log 47r + log r (^1 + y j - — log(?n + 4). 



TTi 

log C(m) ~ log Im ^ — log m as m — 7- oo. 
This proves (1.4). □ 

4.2. Wigner's semicircle law revisited. We can now present the postponed proof of (4.6). The 1- 
point con^elation function Rn{x) can be expressed explicitly in terms of Hermite polynomials, [27, 
Eq. (7.2.32) and §A.9], 

n-l 



n-i 1 „ 

R^{x) = S2i)k{xf + i^n-i{x) / e{x-t)iJn{t)dt + an{x), (4.8) 




where 



and 



M^) = ^e-'^H^ix), H^{x) = i-ire^'f^ie-^'] 

Jo, nG2Z, 

i, X > 



e{x) = <( 0, x = 0, 
X <Q. 



^There are different rescalings of the semicircle measures in tlie literature. Our conventions agree witli those in [27]. 
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From the Christoffel-Darboux formula [42, Eq. (5.5.9)] we deduce 

n— 1 n—1 ^ ^ 

k=0 k=l ' \ )■ 

Using the recurrence formula H'^^ = 2xHn — Hn+i we deduce 

H'^{x)Hn-l{x) - Hn{x)H'^{x) = Hl{x) - Hn-l{x)Hn+l{x) 

and 

kn{x) = -[HI{x) - Hn-~l{x)Hn+l{x) 

2"(n - l)!7r2 ^ 

We set 

r ( ^ K{Vnx) £n{Vnx) - 1 _ 

kn{x) := -= , in{x) := -= , Rn{x) = —=Rn[\/nx) = pn{x) 

\ n Jn Jn 



so that 

Rn{x) = kn{x) +In{x). 

Lemma 4.1. 

lim sup 1 4 (2;) I = 0. (4.9) 



Proof. Using the generating series [-. 2, Eq. (5.5.7)] 

00 , 

we deduce that 

f / Hn{x)dx ) — r = e^^ / e 2 dx = \/27re^^ 



rpn 

"nV-^j— r - e 
n! 
n=0 



n=0 

SO that 



- — — / e 2 H2n{x)dx = — — and / ip2n\x)dx = — ^ ~ const • as n — )• oo. 

(2n)! n! 2"n!7r4 

Using [13, Thm. 6.55] or [ :, Thm. 8.91.3] we deduce that 

svip\ipn{x)\ = 0(n^3^) 

and thus 

sup |an(x)| = 0(?z^T2~4) = 0(?z^3 ) as n — )• oo. 

We set 

^^(x) = / e{x-t)i;n{t)dt. 
Jr 

Using [13, Thm. 6.55 + Eq. (6.26)] we deduce sup^.g]g I F„ (x) I = O ( n- 12 ) . This proves (4.9). □ 
From the above lemma we deduce that 
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Lemma 4.2. 



lim (kn{s)-p{s))wn{s)ds = 0. 




Figure 1. The graph ofkiQ{x), \x\ < 2. 



Proof. Fix c G (0, \/2) so that the interval (— c, c) lies inside the oscillatory regime of Hn{\/nt). We 
have 



\s\<c 



(fcn(s) - p{s))wn{s)ds 



kn{s) - p{s))wn{s)ds + / {kn{s) - p{s))wn{s)ds 



\s\>c 



< sup \knis) - p{s) \ + sup 



\s\<c 



\s\>c 



(fc„(s) - p{s))\ / Wn{s)d, 
J\s\>c 



Using the Plancherel-Rotach formulae ([13, Eq. (6.126)], [37], [42, Thm. 8.22.9]) and arguing as in 
[17, §7.1.6] or [19, §6.1] we deduce that 



On the other hand 



and [42, Thm.8.91.3] implies that 



lim sup \kn{s) — p{s)\ = 0. 



lim / Wnis)ds = 0, 

'^-^°^J\s\>c 



sup I ( kn{s) — p{s) ) I = 0(1) as n — )• oo. 

|s|>c 



□ 



Since Wn{s)ds converges to the (5-measure concentrated at the origin we deduce 

V2 



lim / p{s)wn{s)ds = p(0) 



TT 



This proves (4.6). 
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Appendix A. Gaussian measures and Gaussian random helds 

For the reader's convenience we survey here a few basic facts about Gaussian measures. For more 
details we refer to p]. A Gaussian measure on R is a Borel measure ^rn,a of the form 

ImA^) = dx. 

ay ZTT 

The scalar m is called the mean while a is called the standard deviation. We allow a to be zero in 
which case 

7m.o = <5m = the Dirac measure on M concentrated at m. 
Suppose that V is a finite dimensional vector space. A Gaussian measure on V is a Borel measure 
7 on V such that, for any ^ € , the pushforward C*(7) is a Gaussian measure on M, C*(7) = 

7m(5),a(0- 

One can show that the map 9 ^ i— )• m(^) G M is linear, and thus can be identified with a vector 
G V called the barycenter or expectation of 7 that can be alternatively defined by the equality 
= vd'y{v). Moreover, there exists a nonnegative definite, symmetric bilinear map 

E : X ^ M such that a{Cf = ^(^,0, G . 

The form I] is called the covariance form and can be identified with a linear operator S : V"^ — V 
such that 

where {—,—) : V'^ x V — M denotes the natural bilinear pairing between a vector space and its 
dual. The operator S is called the covariance operator and it is explicitly described by the integral 
formula 

{^,Sr]) = A{^,rj) = / {^,v -m^){r],v -m^)d-f{v). 
Jv 

The Gaussian measure is said to be nondegenerate if XI is nondegenerate, and it is called centered if 
m = 0. A nondegenerate Gaussian measure on V is uniquely determined by its covariance form and 
its barycenter. 

Example A.l. Suppose that U is an n-dimensional Euclidean space with inner product (— , — ). We 
use the inner product to identify U with its dual U'^ . If A : U ^ U is a. symmetric, positive definite 
operator, then 

d^A^) = -— ^^==e-i(^"'"'") \du\ (A.1) 
(27r) 2 Vdet A 

is a centered Gaussian measure on U with covariance form described by the operator A. □ 

If V is a finite dimensional vector space equipped with a Gaussian measure 7 and L : V ^ U is 
a Unear map then the pushforward L^,7 is a Gaussian measure on U with barycenter 

and covariance form 

Si,^ : X ^ E, Si,.^(?7,?7) = S^(Z^r/,L^r/), V?? G ?7^, 

where : [/^ — )• V'^ is the dual (transpose) of the linear map L. Observe that if 7 is nondegenerate 
and L is surjective, then Z*7 is also nondegenerate. 

Suppose (S, fi) is a probability space. A Gaussian random vector on (§, fi) is a (Borel) measurable 
map 

X : § ^ V, V finite dimensional vector space 
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such that X^^fiis a Gaussian measure on V. We will refer to this measure as the associated Gaussian 
measure, we denote it by and we denote by Sx (respectively Sx) its covariance form (respec- 
tively operator), 

^x{ii,i2) = E{{ii,X - E{X)) {i2,X - E{X))). 

Note that the expectation of 7x is precisely the expectation of X. The random vector is called 
nondegenerate, respectively centered, if the Gaussian measure 7x is such. 

Suppose that Xj : § ^ Vi, j = 1, 2, are two centered Gaussian random vectors such that the 
direct sum Xi X2 : S — ® V2 is also a centered Gaussian random vector with associated 
Gaussian measure 

7XieX2 = PXi®X2{xi,X2)\dxidx2\. 

We obtain a bilinear form 

coviXi,X2) ■.V'(xV^2^R, cov{Xi,X2m,^2) = 

called the covariance form. The random vectors Xi and X2 are independent if and only if they are 
uncorrelated, i.e., 

cov{Xi,X2) = ^. 

We can form the random vector E{Xi\X2), the conditional expectation of Xi given X2. If Xi and 
X2 are independent then E{Xi\X2) = E{Xi), while at the other extreme we have E{Xi\Xi) = Xi. 

To find a formula for E{Xi\X2) in general we fix Euclidean metrics {—,—)Vj on V j. We can 
then identify cov{Xi, X2) with a linear operator Cov{Xi, X2) : V2 — Vi, via the equality 

E{{^l,Xi){^2,X2)) = COv{Xi,X2m,^2) 

= {Ci,Cov{x,,X2)cl), vciGvy, 6gv^, 

where ^2 £ ^2 denotes the vector metric dual to ^2- The operator Cov{Xi, X2) is called the 
covariance operator of Xi,X2. 

Lemma A.2 (Regression formula). IfXi and X2 are as above and, additionally, X2 is nondegener- 
ate, then 

E{Xi\X2) = Cov{Xi,X2)S'^\[X2 - E{X2)) + E{Xi). (A.2) 

Proof. We follow the elegant argument in [4, Prop. 1.2]. We seek a linear operator C : V2 — Vi 
such that the random vector Y = Xi — CX2 is independent of X2. If such an operator exists then 

E{Xi\X2) = E{Y\X2) + E{CX2\X2) = E{Y) + CX2 = CX2. 

Since the random vector Xi — CX2 is Gaussian the operator C must satisfy the constraint 

cov{Xi - CX2,X2) = O^^O = Cov{Xi - 0X2^X2) = Cov{Xi,X2) - Cov{C X2,X2) 

To find C we note that 

( ei, Cov{CX2, X2)il ) = E{ (6, CX2)(6, X2) ) 

= i;((c^6,^2)(6,^2)) = Sx,(c^ei,e2) = 

Identifying V2 with via the Euclidean metric ( — , — ) , we can regard 8x2 a linear, symmetric 
nonnegative operator V2 — )• V2, and we deduce that Cov{CX2,X2) = CSx2 = Cov{Xi, X2) 
which shows that 

C = Cov{Xi,X2)S-x\. (A.3) 

□ 
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The conditional probability density of Xi given that X2 = X2 is the function 



P(Xi|X2=a;2)(^l) 



For a measurable function / : Vi — t- M the conditional expectation E{f{Xi)\X2 = X2) is the 
(deterministic) scalar 

E{f {Xi)\X2 = X2)= / f{xi)p(^Xi\X2=a^2)i^i)\dxi\- 

JVi 

Again, if X2 is nondegenerate, then we have the regression formula 

E{f{Xi)\X2 = X2) = E{f{Y + Cx2)) (A.4) 
where Vi is a Gaussian vector with 

E{Y) = E{Xi)-CE{X2), Sy = Sx,-Cov{Xi,X2)S^ICov{X2,Xi), (A.5) 
and C is given by (A.3). 

Let us point out that if X :§—;■[/ is a Gaussian random vector and L : U ^ V is a. linear map, 
then the random vector LX : S — > V is also Gaussian. Moreover 

E{LX) = LE{X), 1:lx{^,0 = ^x{L''^,L''0, G y\ 

where L'^ : V"^ — )■ is the linear map dual to L. Equivalently, Slx = LSxL"^ ■ 

A random field {or fiinction) on a set T is a map ^ : T x (S, ^) — )• M, {t, s) 1— )• £,t{s) such that 

• (S, //) is a probability space, and 

• for any t the function : § — ^ is measurable, i.e., it is a random variable. 

Thus, a random field on T is a family of random variables parameterized by the set T. For 
simplicity we will assume that all these random variables have finite second moments. For any t G T 
we denote by iit^ the expectation of The covariance function or kernel of the field is the function 
Q : T X T ^ M defined by 

CiitiM) = E{{Cti -Mii)(?i2 -Mi2)) = / (61 (s) - lJ-ti){Ct2{s) - 1^-12) dfi{s). 
The field is called Gaussian if for any finite subset F C T the random vector 

is a Gaussian random vector. Almost all the important information concerning a Gaussian random 
field can be extracted from its covariance kernel. 

Here is a simple method of producing Gaussian random fields on a set T. Choose a finite dimen- 
sional space U of real valued functions on T. Once we fix a Gaussian measure d'j on U we obtain 
tautologically a random field 

^iTxU^R, {t,u)^ ^t('ii) = u{t). 
This is a Gaussian field since for any finite subset F C T the random vector 

is Gaussian because the map E is linear and thus the pushforward H*(i7 is a Gaussian measure on 
M^. For more information about random fields we refer to [1, 4, 12, 21]. 

In the conclusion of this section we want to describe a few simple integral formulas. 
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Proposition A. 3. Suppose V is an Euclidean space of dimension N, f : U ^ M. is a locally 
integrable, positively homogeneous function of degree k > 0, and A : U ^ U is a positive definite 
symmetric operator. Denote by B{U) the unit ball ofV centered at the origin, and by S{U) its 
boundary. Then the following hold 



„ / f{u)\dA{u)\ = \ j f{u)\du\ 

-K^ {k + N) J S(U) TT^ Jb{u) 



/(") — —\du\- 
U 7r~ 



(A.6) 



r(i + ^: 

/ f{u)d-f,AH =t^ f f{u)d-fAu) Vt > 0, (A.7) 
Ju Ju 
where d'jj^ is the Gaussian measure defined by (A.l). 



Proof. We have 



fiu)\du\ = ft^+^-^ [ [ f{u)\dA{u)\] = [ f{u)\dA{u)\. 

B{U) JO \JS{U) I K + 1\ Js(u) 



On the other hand 

^ / f{u)e-\^\"\du\ = r t^+^'-^e-'^d^ [ f{u)\dA{u)\ 

ttT Ju vr^ \J0 / Js(U) 

2tt~ \ ^ J Js(U) 2tt~ \ ^ J Jb(U) 

vr^ V ^ / Jb(U) 

This proves (A.6). The equality (A.7) follows by using the change in variables u = t^v. □ 

Appendix B. Gaussian random symmetric matrices 

We want to describe in some detail a 3-parameter family of centered Gaussian measures on §„, 
the vector space of real symmetric m x m matrices, m > 1. 
For any I < i < j define ^ij G so that for any A G 

Cij{A) = Qij = the {i, j)-th entry of the matrix A. 

The collection (6y )i<j<j<m is a basis of the dual space S^. We denote by {Eij)i<i<j the dual basis 
of Sm- More precisely, Eij is the symmetric matrix whose and {j, i) entries are 1 while all the 
other entries are equal to zero. For any A G we have 

The space is equipped with an inner product 

(-,-) :§„xS^^M, {A,B)=tr{AB), yA,Be§m. 
This inner product is invariant with respect to the action of SO(m) on §„. We set 
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The collection {Eij)i<j is a basis of Sm orthonormal with respect to the above inner product. We set 



The collection {iij)i<j the orthonormal basis of dual to {E-ij). The volume density induced by 
this metric is 

\dX\ ■.= J\diij = 2^i^)j\dxi,. 

i<j i<j 

To any numbers a,b,c satisfying the inequalities 

a-b, c, a+ {m - l)b > 0. (B.l) 

we will associate a centered Gaussian measure Ta,b,c on uniquely determined by its covariance 
form 

defined as follows: 

^{^ii,^ii) = a, 'S{Cii,Cjj) = b, yi^j, (B.2a) 

^{^^J,^ij) = c, = 0, Vi<j, k<£, {i,j)^ik,i). (B.2b) 

To see that Sa,fe,c is positive definite if a,b,c satisfy (B.l) we decompose as a direct sum of 
subspaces 

2 j 

With respect to this decomposition, and the coiTcsponding bases of these subspaces the matrix Qa,b,c 
describing 'Sa,b,c with respect to the basis has a direct sum decomposition 

Qa,b,c = Gm{a,b) © Cl(n,y 

where Gm{a, b) is the m x m symmetric matrix whose diagonal entries are equal to a while all the 
off diagonal entries are all equal to b. 

The the spectrum of Gm{a, b) consists of two eigenvalues: (a — b) with multiplicity (m — 1) and 
the simple eigenvalue a — 6 + mb. Indeed, if Cm denotes the in x m matrix with all entries equal to 1, 
then Gm{a, b) = {a — b)lm + bCm- The matrix Cm has rank 1 and a single nonzero eigenvalue equal 
to m with multiplicity 1. This proves that Qa,b,c is positive definite since its spectrum is positive. We 
denote by dTafi^c the centered Gaussian measure on §m with covariance form Xl^ ^ c- 

Since is equipped with an inner product we can identify Sa,fe,c with a symmetric, positive 
definite bilinear form on S^- We would like to compute the matrix Q = Qa,b,c that describes 'Sa,b,c 
with respect to the orthonormal basis {Eij)i<i<j. We have 

Q{Eii, Eii) = Q{iii, lii) = a, Q{Eii, Ejj) = b, Vi / j, 

Q{Eij,Eij) = Q{iij,iij) = 2Q{^ij,S,ij) = 2c, Vi < j, 

Thus 

Qa,b,c = Gm{a,b)e2cl^..y (B.3) 
If I — \a,b,c denotes the Euclidean norm on determined by Sa,fe,c then for 

A = aijEij = an Eii + "V^^^ o-ijEij. 

i<j i i<j 
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we have 



\A\l,b,c = a^al + 2bJ2o'iiO'jj + '^<^Y1 

i i<j i<j 



i \ i / \ j i<j 

= {a-b-2c)'^al + b{tT:Af + 2ctTA^. 

i 

Observe that when 

a-b = 2c (B.4) 

we have 

\A\l^^^ = 6(tr Af + 2ctr A^ (B.5) 

so that the norm | — \a,b,c and the Gaussian measure dVa^^c are 0(m) -invariant. Let us point out that 
the space equipped with the Gaussian measure dr2,o,i is the well known GOE, the Gaussian 
orthogonal ensemble. 

To obtain a more concrete description of Ta.b,c we first identify Xla,b,c with a symmetric operator 
Qa,h,c '■ Sm — ^ Sm- Using (B.3) we deduce that 

Qa,b,c = G{a,h)®2cl^..y 

Observe that 

detQa,b,c = (a - 6)(a + (m - l)&)"'"\2c)(2), (B.6) 

and 

Qaic = Qa',b',c' = G^{a', b') e 2c'l(™-), (B.7) 



where 2c' = ^ and the real numbers a' , b' we, determined from the linear system 

a'-b' = ^ 

a—b 



(B.8) 



a'+{m-l)b' = — r^— nr. 
We then have 

dTa,bAX) = ^(^^.) ^ ^ ^e-^(^-t^'^)2HT) J]dx,,, (B.9) 

(27r) 4 (det Qa^fe^c) 2 i<j 



where 



(B.IO) 



The special case 6 = c > 0, a = 3c is particularly important for our considerations. We denote 
by (— , — )c and respectively dTc the inner product and respectively the Gaussian measure on 
corresponding to the covariance form S3c,c,c- 

If we set Qc ■= Q3c,c,c then we deduce from (B.7) that 

1 

2c (2)' 

where 



Qc^ = Qa',b',c' = G^{a',b') e ^1, 



a'-b' = ^ = 2c' 
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We deduce 



mb' 



1 



1 



m 



b' 



1 



(m + 2)c 2c 2c(m + 2) ' " 2c(m + 2) ' 
Note that the invariance condition (B.4) a' — b' = 2c' is automatically satisfied so that 

(Q^^X,X) = —trX^ r(trX)2. 

' 2c 2c(m + 2)^ ' 



Using (B.6) and (B.9) we deduce 

dT,{X) 



1 



7?i ( m + 1 ) 

(27rc) 4 



. e-^(*^^'"^(*'-^)')2K'?)j]|dx,,-|, (B.ll) 



i<3 



\dX\ 



where 



^„ := 2(2)+(™-i)(m + 2). 
The inner product (— , — )c has the alternate description 

{A,B)c = Ic{A,B):=Ac / [Ax,x){Bx,x)^^\dx\ 

JM.rn, TT 2 

e 2 



= c/ (Aa;,a;)(5a;,a;)- — -^\dx\, \/A,BeSm- 
Jr™ (27r) 2 

To verify (B.13) it suffices to show that 

Ic{Eii,Eii) = 3c, Ic{Eii, Ejj) = c, Ic{Eij,Eij) = 4c, VI < -i < j < m, 

UEij,Eke) = 0, yi<i<j<m, k < £, ^ {k,£). 

To achieve this we need to use the classical identity 

r^Pi . . . T.2pm„-|a;| 



Observe that 



e-|x|2 

{EiiX,x){EjjX,x) — ^|da;| 



vr 2 



vr 2 



m /l\"l-2 



'r(|)r(i), z = ^- ^ 



3, « = J 

1, 



Next, if i < J we have 



vr 2 Jrom - TT 2 



{EijX,x){EijX,x) ^ _rn \dx\ = 4: I xfxj — ^\dx\ = 1. 



Finally, if i < j, k < £ and (i, j) / {k,£), then the quartic polynomial 

{EijX,x){Ekex,x) 

is odd with respect to a reflection Xp 1— )• — Xp for some p = {i,j,k,£} and thus 

{EijX,x){Ekex,x) — ^|da;| =0. 
' vr 2 



(B.12) 



(B.13) 
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Appendix C. Some Gaussian integrals 
The proof of (2.38). We want to find the value of the integral 

4(2^)2 Jsa l<i<j<2 

We first make the change in coordinates 

xii=x + y, X22 = x-y, xu = z. 



Then 



Hence 



(x = V2u) 



detX = x^ -y"^ - z"^, tiX = 2x, ti = lix"^ + y^ + z"^) 

^ [ \x^-y^-z^\e--^'^-'+^y'+^''^V2\dxdydz\ 
Jr3 



2{2n 



(27r)t Jr^ 



1_ I \2u^-y^-z\-^^'+y'+''^\dudydz\. 
7r2 Jr3 



We now make the change in variables y = r cos 6,y = r sin 0, r > G [0, 2tt) and we deduce 



2 



OO / f2lT 



I = — I I ( / \2u^ _ r^le"" dO I rdrdu 



vra JR JO \Jo 

|2n2-r2|e-("'+'''Vdrdu 



Stt 



3 

TT2 Jo Jo 



We now make the change in variables 

TT 

u = tsmip, r = tcosip, t > 0, < (p < —, 

and we conclude 



/ = — J- I / e t (it ) I / 1 3 sin 97 — 1 1 cos ipdip 



(t = y/s, X = sin If) 



4 



1 



e'^'stds 1 ( / |32;2 _ ) = ^ X r(-) X — p = □ 



7r2 \Jo / \J0 / 7r2 

^2/ 3\/3 \/3 

Proposition C.l. Consider the Gaussian measure r„ on $2 with covariance fonn c„. where 



I^--—f detXdTniX) = 2. (C.l) 
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Proof. Using (B.9) we deduce 

dTn{X) = 3 \, _e-K«''-f--n^,^)^ ^ dxij, 

(27r)2(detQa„,fe„,c„)2 i<i<j<2 

where 

detQa„,6„,c„ ^^=^^ 2(a„ - 5„)(a„ + 6„)Cn = St^iSn + '^■tn) =: 5n, 
and, according to (B.8), the parameters a'„, 6'„ are determined by the linear system 



a„ — U 



n a„-b„ 
I I U — 1 



We deduce that 



and since Un — bn = 2cn we have 



bn 



a? - 62 



Hence 

Using the change in variables X := y/UiX we deduce 

Sn 87r2(s„ + 2t„)2 i<r<i<i 



We make the change in coordinates 

xii=x + y, X22 = x-y, xi2 = z. 

Then 

detX = x^ - - z^, tTX = 2x, ti X"^ = 2{x^ + y"^ + z^) . 

2{Sn + 2tn) ^ ^ ' 2{Sn + 2t„) 

^ 4{sn + 2tn){y'^ + z'^) + Unx'^ 

2{Sn + 2tn) 

and 

^ f (an+2t„)(i;^ + z^) + tn^^ 

J„ = 2V2 / (x^ - {y^ + z^))e 2(.„+2*„) dxdydz 

= 2V2(s„ + 2t„)2 / (x^-{y^ + z^))e 2 dxdydz 

{x := x^/2, y := yV2, z := 2\/2) 
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16(s„ + 2t„) 



We now use the identity 



u\-i^^+^'+^^)dudvdw = / -u^e"" x / e^' dt] = vr 



Hence 





poo 


27r 






Jo 










■ X 









3 

7r2 



3 3 

_ , — .. , ^ _ - , _ Sttz (g„ + 2t„)2 

" 1 ^oIj. iOJ.) — 1 ^ 



Using the last equahty in (C.2) we obtain (C. 1) . □ 
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